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A real space diagrammatic method, which is an extension of the Berezinskii technique to problems 
s ! ■ with periodic boundary condition, is formulated to study the density of states (DOS) p(e, (j>) and its 

moments for a one-channel weakly disordered ring threaded by an external magnetic flux <j>. The 
exact result obtained for the average value of the DOS shows that p(e, <j>) oscillates with a period 
of the flux quantum <f)o = —■ However all higher moments of the DOS oscillate with the halved 
period 4p-. The exact expression for the DOS is valid for both weak localization (L 2> I, where L 
is the rings circumference and I is the mean free path) and ballistic (L < I) regimes. In the weak 
localization regime the distribution function of the DOS is calculated, which turns out to be of 
logarithmic normal form. 
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I. INTRODUCTION 



Interference effects in low-dimensional disordered conductors still attract attention from both experimental and 
theoretical physicists, although all main features and a lot of new effects have already been discovered during the last 
twenty years. The prediction of the oscillation in the kinetic coefficients in multiply connected disordered normal metals 
in an external magnetic fieldtl and its experimental observations in a Mg cylinder was a very excellent examination 
of weak localization phenomena, since the coherence of the electron wavefunction during the circulation of a closed 
contour is required to observe an oscillation. The period of oscillation of the magnetoresistance predicted and observed 
firstly was equal to half of the magnetic flux quantum (f>Q — — . Further improvements of the experiments on rings with 
large diameters and small widths gave rise to the observation of a magnetoresistance oscillation with the period (po[ 
0], furthermore in the experiment of Chandrasekhar et ala both periods were observedoO Such a complex magnetic 
field dependence of the magnetoresistance seems to be related with the statistical properties of the sample. 

j» Another development in theory was the prediction of a persistent current in a one-channel disordered isolated 

T — I \ loopfl Due to the analogy between loop and one-dimensional (ID) lattice with a period equal to the circumference of 
■ the loop, a circulating current in rings was suggested which is a periodic function of the enclosed flux with a period 
of 4> . Studies of the effects on the persistent current at finite but small temperature and weak inelastic scattering 
show that both weak-inelastic and elastic scattering do not destroy it.HTLil In experiments, the persistent current 
was also observed. E3~t3 A magnetization measurement was performed in [ |l2| onJV = 10 7 disconnected copper loops 
at T < 1.5K where the electron phase coherence length L v exceeds the loop's circumference L = 2.2/xm, and it 



f 

shows evidence for a flux-periodic persistent current with halved period and 3 x lO -3 ^ 1 amplitude per ring, which 



is remarkably higher than the theoretically expected value for the persistent current per ring, J [ liqjlq jri 



Measurements of the persistent current on single loops (with at lj^st a few channels) in the diffusive^ (L = 8/im 
and I = 70nm, where I is the elastic mean free path) and ballistictfl (L = 8.5/xm and I = Hyum) regimes reveal the 
period 4>q. The amplitude of the harmonic with 4r was measured to be smaller by a factor of 2-3 than that of the 



harmonic in [ 13 



Effects of impurities in most of the theoretical investigations were taken into account by the trsuasfjsj'-aiatrix method 
according to the Landauer expression and by generalizations of this method to n-channel systems JiaclJ'O klihongh the. 
Landauer formula gives full-flux periodicity for all physical parameters, averaging over ensembles of ringsB'EjO'EiTO 
or the calculation of the dynamical current instead of the thermodynamical potentialEjEj were suggested as an 
explanation of the observed halved periodicity. In the process of averaging over different impurity realizations in the 
ensemble, the number of particles in each ring is proposed to be constant, i.e. the persisteet current is assumed to be 
determined by the thermodynamic potential instead of the grand canonical potential£3~ cJ Although there has been 
done a lot of work on the Aharonov-Bohm effect, the existing theories of. nan-interacting electrons stili can explain 
neither the high value of the experimentally observed persistent currentstSLJ nor its diamagnetic sign.113 This can be 
partially connected with the complexity of the experiments, particularly with pdifficulties of the separation of phase 
effects in the rings with relatively large width (larger than the mean free path)li3'Ej from orbital ones. 
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On the other hand, correlation effects may be a reason for the discrepancy between theory and experiment.^ 
Unfortunately, there is still no agreement on the effects of .Co uLom b interaction on the .-amplitude of the persistent 
current. Studies based on spinless electrons in ID continuunJ2SEIK3 and lattice modelscjTa gave controversial results, 
so the amplitude of the persistent current was shown to be increased up to its disorder-free value according to the 
former model, but a Mott-Hubbard metal-insulator transition in the latter model was found to reduce the amplitude, 
or at least a possible increase of the amplitude was negligibly small. There was also the suggestion that correlation can 
change the fundamental period with the magnetic flux and create fractional periodicity in a ID ringES Impurities and 
correlation acting together are again subject of controversy. Thus, considering weak localization corrections in first 
order in the electron-electron interaction to the grand canonical potential,c!J a persistent current with a period of ^ 
and an amplitude of ~ ^jir, corresponding to the experiment,£3 was obtained, while Monte Carlo simulation on a ID 
Luttinger liquidEj resulted in a persistent current with period </>o and with an amplitude decreased through interaction. 
Resumming the existing results it can be said that neither the noninteracting electron model nor models of correlated 
electrons yet gave satisfactory answers to the questions put forward by the experiments. These concern the period 
(under what condition both periods or only the halved period are observed) , the amplitude of the persistent current 
and its diamagnetic sign (for the correlated electron model the diamagnetic sign requires an attractive interaction 
between the electrons). 

To prevent the interference of the orbital effects in the presence of an external magnetic field with the phase effects, 
the width of the ring should be chosen as narrow as possible, i.e. a one-channel ring with random impurities seems 
to be an ideal tool to study the Aharonov-Bohm effect. However, interference effects in ID disordered systems as a 
result of the coherent backscattering processes are stronger irrespective of the degree of randomness. The diffusion 
approximation, which was used in previous studies, is not acceptable in f D systems even for the case of weak disorder. 

Approaching the problem thoroughly, we use in this paper a weakly disordered noninteracting electron model-and 
construct for it a new exactly solvable diagrammatic method, which is an extension of Berezinskii's methodcilea to 
the problem with periodic boundary conditions. Within this model, we sum up all impurity scattering diagrams in 
the framework of the Born approximation. 

In Sec.|l| we describe the method. We calculate the density of electronic states (DOS) in Sec. III. Indeed an average 



value is not enough to describe the observable parameters in low dimensional mesoscopic systems. It is well known 
that the physical parameters of a mesoscopic systein-with size L satisfying the condition / < L <C L v fluctuate from 
sample to sample, i.e. self-averaging is violated. c^TEI At T-^Ji all sufficiently large systems become mesoscopic. In 
this case, high moments give a considerable contribution,c3~tj which results in strong differences between average 
value and typical one of the observable parameter, i.e. the average value loses its significance to characterize the 
experimental observation. In Sec.^ the diagrammatical method is applied to find the feth moments of the DOS, 
(p k (e,(f>)). The obtained equations for (p k (e,<fi)) show that, in contrast to the average value of the DOS, all higher 
moments oscillate with the halved period, Although the structure of the equations is complicated, the latter can 
be solved for the weak localization regime when the condition I <C L is satisfied. This procedure is given in Sec]v|. 

The zeroth (not depending on <fi) and the first (oscillating with 4r) harmonic of (p fc (e, (f>)) are studied explicitely in 
this section. Both are shown to increase with k as exp(/c 2 ), which gives rise to a logarithmic normal distribution. 



In section VI we conclude our results and discuss possibilities to extend our approach to related problems. 



II. DESCRIPTION OF THE METHOD 

We consider here a one-channel disordered ring with length L = 2irr, threaded by a magnetic flux <p through the 
opening. The Hamiltonian of the system can be written as 

where 4>q = — is the fundamental period of a flux quantum and m* is the effective mass of an electron. The potential 
Vi m p of randomly distributed impurities is considered to be weak, so that scattering processes can be studied in the 
framework of the Born approximation. Below, we use the spatial variable x = np instead of the angle tp. 

Our aim in this section is to construct a diagrammatical method for the calculation of the average values of the DOS, 
(p(e, </>• L)), and of its moments (p™(e, 4>\ L)). The bracket (. . .) denotes the average over the impurity configurations. 
Expressing the DOS by means of retarded (Gr) and advanced (Ga) Green's functions (GF) as 

p{e,<f>;x) = --ImG R {e,cf);x,x) = — [Gvt(e, 4>;x,x) - G R (e,<f>;x,x)] (2) 

TT Z1TI 
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the n-th moment of p(e, <fi; x) can be given by 

(p k (e,<f>;x)) = (^rE (^j(-^) l (G l R (e,^x,x)G k A - l (e, ( t ) ;x > x)) (3) 

The Berezinskii diagram techniqueE3@ is applied to calculate the average value of a single GF and higher-order 
correlators (G l R G A ). In contrast to strictly ID disordered wires, quantum corrections to the DOS of a ring turn out 
to exist even for the weakly disordered limit due to periodicity. 

As for infinite systems, we consider as starting point a free particle with wave function ip p (x) oc exp(zpx), where 

the momentum can assume arbitrary values, leading to a continuous spectrum e p = 2^v(p — ) 2 - The "bare" GF 
G° R A can be calculated easily: 

G RA e^x,x')=( f*^^*^*^^^-^^ (4) 
J 2tt e — e p ± irj v(e)n 

where L is the circumference of the ring and the parameter 77 is introduced phenomenologically to model inelastic 
processes, which result in a blurring of the energy levels. v(e) — yj ^ and hp(e) — y / 2m*e are the velocity and the 
momentum of an electron with energy e, respectively. Notice that Zeeman splitting has not been taken into account 
here. For an electron with spin s = ±i, it would throughout the paper lead to a shift of the energy as e — ► e — sgp R B 
(where g is the gyromagnetic ratio of the electron and ps is the Bohr magneton) . 

This bare GF, however, is not yet the real GF for an electron in a ring without impurities, since it does not reflect 
the finite size of the system and the periodic boundary conditions. These are taken into account by allowing the 
particles to make arbitrary revolutions around the ring, which leads to the expected quantization effect. According 
to this prescription, the GF for a clean ring G R A is 

G° RtA (e, 0; x, x') = G% A (e, cj>; x, x') + G% A {e, 0; x, x' + L) + G% A (e, <f>- x, x' - L) + G% A (e, <f>; x, x' + 2L)... (5) 

One may verify this approach by calculating the DOS of a clean ring in a mag netic field po(e» from Eqs.(§), (@), 
and (|): 



p (e, <j)) = po + 2p cos(p(e)Ln) cos(27r— n)e "<"«) Ln (6) 

where po — l/(Trhv(e)) is the DOS for a pure and infinite ID system. As 77 — > 0, Eq.(^|) displays the discrete behavior 
for the DOS of a clean ring. 

In the diagrammatical technique, the retarded (advanced) GF of an electron moving in the field of randomly 
distributed impurities is represented by an ordinary (double) continuous line in real space, which goes from point 
x to x' after multiple scattering on a given impurity configuration, realized by the potentials V{xt) with impurities 
placed at the points {xi}. For the DOS and its moments it suffices to consider x = x'. Since the "bare" GFs [Eq.(@)] 
between two successive scatterings have factorable structure, the coordinate dependence can be transfered from the 
lines to the vertices. Averaging over the random Gaussian potential leads to a pairing of the impurity vertices. Their 
strength is measured by the inverse forward (backward) scattering length, 

1 2 f°° 

U(x)dx and = — ^ / U(x)dx cos (p(e)x) (7) 



l+(e) h 2 v 2 (e) Jo l~(e) h 2 v 2 {e) 

For the Born approximation to be applicable, the correlator U(x — x') = (V(x)V(x')) should have a width much 
smaller than the mean distance between impurities -, and - <C ^ ■ In the extreme case of a white noise potential 
U(x — x') — ell 2 J^nL-oo ~ x ' + n ^ J \ tli e two scattering lengths become equal, l~ = l + = 21. Given in Fig|l] are 
the essential vertices selected according to the condition ppl 3> 1, with pp and I being the Fermi momentum and the 
mean free path, respectively. Although the "bare" GFs depend on the direction due to the magnetic field, the internal 
vertices in Fig.[y do not differ from those of Berezinskii. All dependence on the magnetic field is transfered from the 
lines to the external vertices, which are shown in Figj^. 





a(a') b(b') c(c') 
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FIG. 1. Contributing internal vertices. Vertices a',b', and c' have the same form as a,b, and c, but double lines instead of 

TO- 



single ones. The following factors correspond to the vertices: a, a' ~ (^4- — jjr); b,b' ~ c,c' ~ j^; d ~ ip-; e ~ ±r exp(-i^l) 
f ~ pr exp(- 




FIG. 2. The external vertices. One incoming vertex and one outgoing vertex are attached to each continuous fermion line 
characterizing one Green's function. 



As an example, a simple diagram contributing to (Gn(e,cf>;x,x)) is drawn in Fig.^a. For convenience, we cut the 
diagram at point x and straighten the lines, which results in Fig.||b. Each diagram for a single GF is then characterized 
by m pairs of lines returning to x and n through-going lines. Since the bulk of each diagram (i.e. after the removal of 
the external vertices) does not depend on the direction, n is the sum of right-going n + and left-going n~ lines. For 



correlators (G R (e, </>; x, x)G^~ (e, <j)\ x, x)), we have to distinguish the number m of returning lines on the l.h.s. and 
the number m! of returning lines on the r.h.s., and in addition to introduce m, m', and n for the advanced GF. An 
example for this case is shown in Fig.^Jc. 

In contrast to Berezinskii's technique for an infinite ID system, the external vertices depend on the direction. Also, 



the diagrams carry a factor exp[27ri-J-(n + + ; 



")] from the through-going lines. The number of pairs on the 



two sides of the cutting line, which we denote by m and ml for the retarded GF and by m and m' for the advanced GF, 
may in general differ by ±1 (i.e. m — m' = —1,0, 1 andm-m' = —1,0, 1) as in Sec. IV. However, for (GR(e,(j>;x,x')) 
in the regime of weak disorder as in Sec. HI, we have always m — m! . 
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a) 



b) 



c) 



FIG. 3. A simple diagram contributing to the single GF (Gr) a) resembling the motion of the electron around the ring and 
b) after cutting at point x. The numbers characterizing this diagram are m — 1 and n = 1. c) shows a contribution to (GrGa)- 
Here, m = 2, m' = 1, n = 0, m = 1, fn! = 0, and n = 1. 
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a) 



b) 



c) 



d) 



FIG. 4. The different possibilities to attach external vertices. The cases b) and d) contain oscillating factors exp(±2ip(e)x), 
which is connected to the fact that the numbers of line pairs on the left and right hand side differ by ±1. The cases a) and c) 
have m = m' and hence do not contain such factors. 



III. THE DENSITY OF STATES 



The diagrams for the DOS do not exhibit the full complexity presented in the previous section, since they contain 
only retap. ed GFs. Consequently, we can omit vertices d),e), and f) of Fig|l] and set m = m! . Following Berezinskii's 
method J_j we denote the sum of all diagrams with m pairs of returning lines and n = n + + n~ through-going lines 
by Qo(m, n;x — x' — L), where n + (n~) is the number of right-(left-)going lines that cross the whole diagram. Such a 
diagram is shown in Fig.^i and [|b. Since the magnetic field dependence has been extracted from Qq, it only depends 
on the total number of through-going lines. The condition m — ml restricts the possibilities to attach the external 
vertices to the cases a) and c) of Fig||. The average value of the retarded GF can be expressed in terms of the kernel 
Qo as 



«^*-»~j55E £ e [("I") ( B ; l T) + ( m r)("'; , .r) -*»•*•* a-* 



exp{zp(e)i(n + + n ) ^—(n + +n ) — 2iri-^-(n + — n )}Qa(m,n + +n ; L) 



(8) 



The two products of binomials in the bracket of Eq. (|^) characterize the different possibilities to insert the n + right- 
going and the n~ left-going lines between m pairs and correspond to the cases a) and c) of Fig.|| respectively. For 
m = n + = -nT = 0, these two possibilities are degenerated into a point-like diagram. To avoid double counting in this 
case, the third term in the brackets has been added. The combinatorial factor, corresponding e.g. to the configuration 
in Fig.[|a, can be obtained as follows: n + lines can be distributed at m + 1 positions, before each of the loops on the 
l.h.s. and after the last loop, whereas the n~ left-going lines can be inserted at m positions before each of the loops on 
the r.h.s.. Denoting the number of lines at a given position with n i , we have the restrictions n + = +n\ ■ ■ ■ + n^ n+1 
and n~ = n± + ■ ■ ■ + n~ n . 

Summing over all these possibilities gives 



E 



s. 



+ 



{«+} {n-} 

The expression for the DOS can be obtained by combining Eqs.(0) and (ph: 



n- ,n t +n 2 +-+n m I M m -\ j [) 
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m=0 n=0 k=0 



(10) 
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cos(p(e)Ln) exp( r-rLn) cos(27r — (2k — n))Q (m, n; L) 

v(e) 0o 



5 



To find QQ(m,n;x — x'), we shift x infinitesimally and examine the different impurity vertices that pass through 
the point x. The result is a differential equation for Qq, 



d r (2m|n) 2 n m(m + n)-, . 1 2 , 

— g (m,n;x) = -[ — h — H p jg (m,n;x) - — m Qo{m- l,n + 2;x) 



ill) 



where the vertices a),b), and c) in Fig.ffl contribute. For a) the number of possibilities to be included is 2m + n, for b) 
^(2m + n)(2m + n — 1) and to include cj there are m(m + n— 1) ways that do not change m and n and m 2 possibilities 
that decrease m by 1 pair and increase n by 2. The process of construction of this equation is illustrated in Fig.|| 



+ 2 



+ 



+. 



FIG. 5. How to include the internal vertices to Qq. Shown on the r.h.s. is the inclusion of vertex b) of Fig.jl] and the inclusion 
of vertex c) without and with changing of m and n. 



The kernel Qo satisfies the boundary condition 

Qo(m, n; x — x' = L = 0) = 5 m q 



(12) 



which expresses the absence of scattering for a ring with infinitesimal small circumference. The Eq.(|lT|) for Qo(m, n; x) 
can be solved exactly as shown in Appendix |a|: 



Q (m,n;L) =exp{-— p(2m + n) 2 - —m(m +n) - ^p"} 



m\ m\(j+n— 21! ,„ . r L , „ .-. 

1 J- (2j + n-l exp{-jj + n-l } 



(13) 



Equations (|10|) and ( |13[ ) give a complete description of the DOS for a one-channel ring in an external magnetic field. 
They are exact in the regime of weak disorder, when the condition pfI^ 3> 1 or epr 3> 1 is satisfied (pf and tp are 
Fermi momentum and Fermi energy, respectively). Further, since the circumference L of the ring may vary within a 
large range [L < P 1 to I <C L < oo), the results cover both the weak localization and the ballistic regimes. For the 
weak localization regime, when L ^> max(/ + ,^ _ ) the amplitude of Qg(m, n; L) decreases rap idly with m and n due 
to the exponential prefactor. Keeping harmonics up to n = 2 in Eq.po|) and using Eq.( [A6| ), we obtain the leading 
behavior for p(e,4>) 



. r 2L 2L L . -| . i L ,r . . . . , e> . 

p(e, 0) =p [1 - — CX P(-— - p )J + 2 Po ex P( _ ^+ ~ 21=) |cos(p(e)Z) cos(2vr— )e 

L 



+ e iF 



cos(2p(e)i) cos(47r— )e 

00 



(14) 



We see that the main quantum correction to the average value of the DOS oscillates with a period of 4>q. The 
amplitude of this contribution decreases exponentially with the impurity strength or with increasing L, so that p = pa 
for L — > oo. 

The ballistic regime is realized for L < min(Z + ,^). Keeping terms up to first order in m, the DOS can be 
approximated in this limit by 



L N+ <£> L N ~ <j) 

p(e, <f>) = Po( e 7 4>) — Pott / " 2 cos(p(e)Ln) cos(27r— n) — po— } n cos(p(e)Ln) cos(2-7r— n) 



2p p E( fc + C0S (p( e ) Ln ) cos(27r-^-(2fc - n)) 



N~ n 



(15) 



G 



where N + [y 2£t] and AT « [^-] , and the DOS of a clean ring po(e, <jj) is given by Eq.(^). In the ballistic regime, 

L is of the same order of magnitude as , hence may be rather small integers. Therefore the oscillation with the 
full flux quantum <pQ dominates in the ballistic regime. 

In the absence of backward scattering (l~ = oo), Qo(m, n; L) is greatly simplified: 

Qo(m, n;L) = exp{~-^n 2 }8 m . (16) 

Substituting Eq.@ in © and using lim m ^o ( m+ m-i = ^n,ki we can express the DOS as 



cxp(-ip(e)L - i2ir(f>/ </> + iz + ^feL) - 



1 



: — — — , l , , — — — : — — ~ 7 + < •'•} 



(17) 



exp(-zp(e)L + i2-K(f>/tpo + iz+ ^fcjL) - 1 

From Eq.(^) one sees that forward scattering coherently shifts all energy levels. The value of this shifting is random 
with Gaussian distribution, with a typical value of 7V yx") where r + is the relaxation time due to forward scattering. 

The level repulsionS in ID disordered systems therefore is only due to backward scattering. The width of this level 
broadening of the averaged system is much smaller than Dingle broadening for the weak localization regime, whereas 
the two mechanisms can have comparable effects in the ballistic regime. 



IV. HIGHER MOMENTS OF THE DOS AND DISTRIBUTION FUNCTIONS 



Exact calculations show that the average value of the DOS oscillates with a period of the flux quantum 4>q. To 
understand .the reason for the experimentally observed oscillation of the persistent current in a sufficiently large ring 
[L Z)E3ll3 with the halved period, we calculate here higher moments of the distribution of the DOS. According to 
Eq.(||), we have to determine the correlators (G l R G k ^ 1 ) for the fcth moment. In contrast to the Berezinskii technique 
for strictly ID systems, the correlators here are characterized by only one block Q. An example is shown in Fig.^c. 
Each diagram contributing to {G l R G k ^ 1 ) consists of I retarded and k—l advanced lines. For the zth retarded (advanced) 
GF, we count the number of left loops m,i(jf%i), of right loops to^(to^), and of left- and right-going traversing lines 
m = n\ + i\ [rLi = 7i+ i +n~i). (Notice that the index i of ra here denotes the number of the GF, whereas it was used 
for the position within one GF in the previous section.) For the different fermion lines, we now can attach the external 
vertices in four different ways, as shown in Fig.|j. This allows m, L — m! i = —1, 0, 1 (m, — lfi! i — — 1, 0, 1) for the each 
retarded (advanced) GF. It turns out that the block Q depends only on the total numbers m = mi + rri2 + ■ ■ ■ + mi, 
to = toi + TO2 + ■ • ■ + TOfc-i and similarly m' ,rff ,n,n[ |5l| . It is clear from these considerations that the difference 
between ro_and to' is always restricted by — I < to — m' < I. This is automatically taken into account by the mixing 
coefficiently tpi(m, to' ; n + , n~) . This coefficient is the generalization of the term in brackets of Eq.(||). It counts the 
different possibilities to attach external vertices (see FigJJ) and to distribute the through-going lines between the 
loops. Using Eq.(|^), we can write it as 

(Pl(m,m';n + ,n~) = ^ 5 m . mi+nl2 + ... +mi 2J S m ',m' 1 +m! t +-+m> l ^2 S n+,n+ +n++-+n+ 5 n-,n-+n- + -+n-[ 

{mi} {m^} {„+} {„-} 
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mi ,m . 



[ (m l + nj\fm i -l + n i \ fm l -l + nj\fm l + n i \ , "i 

[( ui i )( mt-1 ) + [ mi -l ){ to, J-W^A-o]} 



(18) 



Separating the exponential factors in Eq.jlq), we can write ipi as 

/ / 4- — \ — 2ix(m' — m)(p(e)+i— V^r) ~ i i + — \ /in\ 

ipi(m,m;n . , n ) = e v ww ' (pi(m,m ;n ,n ) (19) 
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with 



m,vai+m 2 A hni; / y °m' ,m' 1 +m' 2 -\ Ym' / y ° n + n +. 

i-m'i} {«+} 

I 



/ j n ,n x 

K") 



+n„ H hn. 



rrii ,m' +1 



n{- 



m,i — 1 + 77^~\ /m.j — 1 + Ji; 



,m' — 1 



rrij + 77^ 



///, I / V ///, 1 / \ III; 

nii + nj \ ( mi — 1 + JiT" \ / TO i — 1 + nt \ (fn%-\- fh 



m l + n i 
m. 



(20) 



m, - 1 



rrii - 1 



rn. 



•*nn,0* n + *n-,o]} 



The mixing factor for the advanced GF is obtained from complex conjugation of (pk—i- 
Now we can express Eq. (0) in terms of Q and the mixing factors 



oc oo oo oo oo oo oc oo 



(/M;^)>^)1EEE E E E E E V/ 



1=0 m=0m=0 m'=0m'=0 n+=0 n~=O^T = o^ = o 

e 



i2TT-^(n + -n +n+-n-) e ip( e )L(n++n--n+-n-) e -^ijj(n + +n-+n++n-) 



-2ip(e)x(m' —m—m! +m) ^jjj) ( m ' —m+rnf -rn) 



(21) 



<£>l(m, m ; n , n 



, + , , — r — > n /m, rn', 77+ + n 

•" " )tpk-i(m,m ;n+,n )Q( _, — r — 



Here, the first three exponential factors come from the external vertices and from the revolutions around the ring. 
The last two exponential factors were separated from the mixing coefficients [Eq. ([!{])]. Below, we shall see that the 
last exponential factor is canceled by a contribution from Q. 



A. Equation for the central block Q 



As noted above, the central block Qi—' _ 

\m, m , n 



x ) is defined as the sum of all diagrams with m and m pairs of 



returning lines on the left side, m! and m! pairs on the right side and n and n through-going lines, coming from 
retarded and advanced GFs, respectively. An equation determining Q can be constructed according to Berezinskii's 
idea by attaching all possible vertices, given in Fig.[j], to the existing block, while avoiding the formation of unconnected 
electron loops. Careful analysis of all these possibilities gives the equation 



d 

d.r 



(m— 1, m' — l, n + 2 
m , n 



(m, m! 
\m, to , 

— —mm U\ 

l~ V 777, 

1 4a»„/m — 1, ml, n \ 

+ — mme^ Q _ -, _,' _ x) 
\— \m — 1, m , n J 

1 / . w_ . _n „ /777 + 1, 777,', 77 

+ — (m + n)(m + 77)e "<«> Q _ . , _/ _ 

\— v 7 \777 + 1, 777 , 77 

+ -^(777 + 77)777'e-^TQf m .I l ' 
l- y ' V 777, 



——(2777 + 77 — 2m — n) 2 + —777(777 + 77) + —777(777 + 77) H — T^C 7 " 1 + n) q(. 



m, rn , 77 
777, rn! , n 



1 , _ m. 777 , 77 

^\777 — 1, 777 — 1, 77 +2 



+ ^-m'(m + n)e ^q(. 

I- ' \r 



1, 77 



7/7 
■/ 



77 

1, 77+2 



1 



, , „/777, 777 

me " <e) y , 

Vm, m 



1, 77 + 2 
1, 77 + 2 



The block Q is subjected to a similar boundary condition as Qq in the previous section 



(22) 



fm, m , 77 
vm, m', 77 



(23) 



which states that for an infinitesimal ring there can be no scattering. The first coefficient on the right hand side 
of Eq.([2^) contains contributions from the vertices a),a'),b),b'),c),c'), and d) in Fig.0. Vertices a),b), and d) can be 
attached in (2m + n),| (2m + 77) (2m + 77 — 1), and (2m + 77) (2m + 77) ways, respectively, the coefficients for a') and b') 
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are the same as for a) and b), with the replacement {to, n} — > {m, n}. For vertex c) we have again to distinguish two 
possibilities as in Sec.|l|. We have to(to + n — 1) ways to attach it without changing to, to', and n; and toto' different 
ways with changing {to, m'n} — > {to — 1, to' — 1, n + 2}. The latter kind of insertion of the vertex c) in Fig.[j] and its 
counterpart for the advanced GF give the second and the third term on the right hand side of Eq. ( f22[) . The inclusion 
of vertex e) reduces to and to by 1. The insertion of vertex f), however, can be done in four different ways which are 
shown schematically in the last four blocks of Fig|| 



I 



3 
+ i 



+. 



FIG. 6. Different possibilities to attach vertices to the block Q, as the left side is shifted. The first plot shows the inclusion 
of vertex c) of Fig.jl] with decreasing m and m' by 1 and increasin g n by 2. The other four blocks show the different possibilities 
to include vertex f) and correspond to the last four terms in Eq.(E2|). 



Trying to solve Eq.(|22|), one may begin by substituting Q = m i{ m i\ Q and then introduce new variables M — 2m + n, 

M = 2m+n, M' — 2m' +n, M = 2m' + n. As a consequence, M' and M appear as fixed parameters in the differential 
equation for Q. But still then, Q depends on the five variables M,n,M,n, and x. Under these circumstances, looking 
for the general analytic solution of Eq.(|2^), one meets with enormous difficulties. Before studying an asymptotic 
approximation of the problem, we make some simplifications of Eqs.(21) and (j22|). The exponential factors in Eq.([2T 
can be removed by the following substitution: 



/to, n N { 2 
^ \to, to , n ) r 1 v 



2r\x 
to 



(to + to — to' 



_ _/to, to', n 
' ' \to, to , n 



(24) 



The equation for Q has the same structure as Eq.([22|), only the exponential factors are dropped and the first term 
on the r.h.s. of Eq.(]22|) acquires another contribution — ^y(to + to — to' — to') to the prefactor of Q. 



From the structure of the internal vertices in Fig.[l] one sees that the condition 

(m — to) — (to' — to) = 



(25) 
) A 



is satisfied for arbitrary cross sections. (The same condition also applies to the strictly ID problem, see [ |5l) 
corresponding symmetry of Eq. ( p2[ ) confirms this condition. 

In the regime of weak disorder, ppl ^> 1, Eq.(pl|) for the fcth moment of the DOS contains terms that strongly 
oscillate with the particle energy (n + + n~ ^ n + + n~) apart from smooth ones (n + + rT = n + + n~). To neglect 
the strongly oscillating terms, we choose only those terms in Eq.(|l]) that satisfy the condition 



n + n = n + + n 



Eq.(Elf) now is simplified to 



oo oo oo oo 



m— m— m'— n— 



where 

$^ (to, to, to' 



m', 



r*. It II / 7 \ 

j ,m ,n,n) = 2^ l^[i) e 'ipi(m,m; 

1=0 n+=0^+ =0 



(26) 



to + to', n 



x^j $fc(m, to, to', to — to + to', n, n) (27) 



n + , n - 



n + )(p k -i{m,m ;n+,n - n+) (28) 



Expressions (27) and ( p8[ ) show that the dominating contribution to (p k ) does not strongly oscillate with the energy. 
Unlike the averaged DOS, the first harmonic of all moments oscillates with the halved magnetic flux 
In the following section, we solve Eqs.(g3)-(p8h for large rings, with L ;» max{i + , 
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V. DISTRIBUTION FUNCTIONS FOR THE WEAK LOCALIZATION LIMIT 



The equation for Q is simplified considerably in the limit of large rings, L ^> max{( + ,r}. For this case, we can 
assume that the electrons are quasi-localized and that the wave function overlaps around the ring are small, similar 
to a tight -binding model. Diagrammatically, this means that the electron loops emerging from the l.h.s. and the r.h.s. 
of the diagram almost nevet-rjeach each other, since they have a characteristic size of £ <C L. (The localization length 
for an infinite ID system isElo £i£> ~ 4l~ .) 

As a consequence, we can for large rings neglect those inclusions of the vertices c),c'),e), and f) that directly connect 
the loops on the r.h.s. with those on the l.h.s.. Corresponding to this is the neglect of the terms 2,3, and 6-10 on the 
r.h.s of Eq.(]25). Now, Q can be factored as 



'TO, 



m 



x J = Q*(m, m, n; x)Q*(m! , to — to + to' ,n;x) 



\to, to — m + TO , n 
where the factors are defined through 

dQ* (to, to, n\ x) r 2?7 2 2 1 1 1 -, 

= — — -r(m + to) + — (to — to) + — m(rn + n) + — rnhn + n) + — n \Q (to, to, n; x) 

ax v ( e ) » » ' i 

+ — totoQ*(to — 1, fn — 1, n; x) + ~( m + n)(m + n)Q*(m + 1, m + l,n; x) 



(29) 



(30) 



Apart from this simplification, the limit £ ^> max(( + , implies to, to, to' n. 

Note that Eq.(|3(]) can also be obtained from Eqs.(p2|) and ( |2^ ) by neglecting to' and to'. The non-entanglement 
mentioned above has a second consequence: The remaining contributions change to and to simultaneously by ±1 
(due to vertices e) and f) of Fig.|l|), or conserve both to and to, as in Berezinskii's approach to strictly ID systems.EJ 
Therefore we can adopt to = to, which further simplifies Eq. (J3Q) : 



l~ - — - — = — \4rjT~m + n + 2to(to + n)l Q*(m, n; x) + m 2 Q*(m — 1, n; x) 

dx 



(31) 



+ (to, + n) 2 Q*(m + 1, n; x) 

Here, r~ is the inelastic scattering time with respect to backward scattering. The boundary condition for Eq.(]3l|) is 

Q*{m,n;x = 0) = 6 mfi (32) 

Th e co mbined mixing function given in Eq.(^8f) is simplified for large to, to, and to' in Appendix^. Substituting 
Eqs.(B6) and (E9) into (pTn, we get a comparatively simple expression for (p k ): 



OO OO OO k / 



2 

m=0m'=0n=0 Z=0 
22n m /c+2n-2 



2/ 



I J \m — m' J \m — m 1 



2k- 21 



(33) 



r(z)r(fc - i) 



(to, n; L)Q*(m , n; X) 



As we emphasized, the diagrammatical structure of the block Q* demands its dependence on one parameter to instead 
of two (m and to). Hereb y th e sum ove r to is removed. The summations over I and to' can be done as described in 
Appendix g. Using Eqs. © and @, we get 



OO OO 



i hi a T w k 2l ~ k (k- l)T 2 (2fc) ^ ^ 



22n e -^-™ TO fe+2ri-2 



fc(2fc-i) r 5 (fc) ^ ^ r 2 (n + i) 

(p k (e, 0; L)) n=0 + (/> fe (e, 05 £))n=i cos 2 (^) + 



[cos(27r- 



1 2n 



(to, n; L) 



(34) 



Eq.(|3l|) for Q*(m,n;x) was solved approximately in [ |32| for arbitrary n. Here we shall study this equation for the 
zeroth and first harmonics [n = and n = 1) in detail. 
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A. Zeroth harmonic contribution to the DOS moments 



The zeroth harmonic of (p k (e, 4>; L)) in Eq.@ contains Q*( ro, n — 0;L). Laplace transforming Eq. (pl|) , written 
for n = 0, with respect to x and using the boundary condition (p2[), we get 

(35) 



(A + Sl m)Q* (m; A) - <5 m , = m 2 [Q* (m + 1; A) + Q* (m - 1; A) - 2Q* (m; A)] 



Here, s% = Aj]t and A is the parameter of the Laplace transform. This is an equation for the right hand side in the 
Berezinskii technique with an open boundary condition and it was solved in [ |53| , pT[ |. Here, we give only the result for 
Q5(si,m; x): 



Ql{s u m;x) = 2(msi)2Ki(2(msi)- 



2(msi) • 
2ni 



dA^i 



i-X T 2 {-iX) 



4-fe- (A +1) ^K_ lA (2(m Sl )3) (36) 



After substitution of this solution into (/3 fe (e, <j>] L)) n=0 in Eq.(^J), the summation over m can be transformed into an 
integration, which is done easier. Some mathematics results in the following form of (p k (e,cf>;L)} n= o: 



(/ M;L )) n=0= (>£) 



PQ . k 2(*-l)r(2A) ,_ k f k 



k(2k-l)T 5 (k)' 



dA e -4^(^n) 2 -^7 2 r 3 (i 



2?ri 



£ — A 

/\|2 



r 2 (-iA) 



w 2fc + l + »A w 2fc-l + a ^ 2 



2 r 00 dz' _^ z -2 

- / «-* r(fc-w') 



(37) 



dz e ar 



i3 /* 1-iz-iz' s \p3 ^ 1-iz+iz' 



«z + iz ) > 



, 27ri (z + z' — z)(z — z' — i) r 2 (— iz — iz')T 2 {—i 

where 7 = ^- In ^- > 0. The second term in the bracket of Eq.(|37j) has a saddle point at Ao = i"f and simple poles 
at the upper half-plane: Ai = i, A2 = i(2k — 1), i(2k + 1), . . . . The integral over z in the third term in the bracket 
contains again the saddle point at zo = 17 and poles at Z\ — ±z' + i, Z2 = ifc, i(fc + 1), . . . . For 7 < 1, the main 
contribution to both integrals is given by the saddle points. As a result we get 



(p fe (e,0;L)) n=o =^(2 Sl ) 



1 _ fc T(2fc-l) 
F(fc) 



(38) 



Such a result has been obtained for the infinite ID disordered system.il Transforming the semi-invariants in Eq.(|3 
to moments and using the inverse Mcllin transformation 



W{ P ) = 



1 

2nd 



a+zoo 



/a— zoo P 

the following inverse-Gaussian distribution function is obtained: 



dfc 

~~k+ 



W n =o(pJ = ( - 3 ) exp(-2 — yyr 



irp° 



PPo 



(39) 



(40) 



For 2r\T 3> 1, the most probable or typical value of p is equal to po, whereas for 2r\T 1 it shifts to lower values 
and becomes equal to p typ — 4n ^ po ■ 

When 7 assumes intermediate values, i.e. 1 < 7 < k, the essential contribution to (p k ) n= o comes from the saddle 
points of the third term in the bracket of Eq. (H?)) and the contributions from the poles of this term cancel the other 
term, resulting in 



(p k (e,frL)) n=0 



p k T(2k - 1) 



,r k i-e-iM 1+ ^>r«(i£n 



2 fc - 1 r(fc- i)r(fc + 1) 7rL(i- 7 ) 2 r 4 ( 7 ) 



r(fe + 7 )r(fc-7) 



(41) 



This expression shows that high moments of the DOS for intermediate values of 7 increase with k; however the increase 
is not so rapid, Eq.(41) has an additional factor ^ compared to Eq.(38). 
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For 7 satisfying the condition 7 > k, the leading contribution is given by the pole zi = ik and 

,k 



P0\k 



(pX(e,c P ;L)) n=0 = (^) 



l- 



2vLk(k- l)(2fe-l) T 7 (k) 



r 2 (2fc) r6( . 



1, 



e2i 



(42) 



The last expression is valid for arbitrary small values of the dissipation parameter (77 — > or 7 — > 00) with r\ <C 
-^rr exp(^p^). Eq.([i"2"|) shows that the zeroth harmonic of the fcth moment of the DOS grows with k as exp(fc 2 ). Such 
rapid increasing of high moments of {p k ) n =o has been firstly obtained by WegneiCj and it is a characteristic feature 
of the logarithmic normal distribution of (p k ) n =o- The distribution function for the zeroth harmonic term can be 
obtained using Eq. (|39|) . For large values of the DOS, satisfying the condition p > ^-exp(p-), the dominating saddle 
point yields again a logarithmic normal distribution: 

p^2i" i„ 2p r,W2Z 



W n=Q {p) 



_ln^ 

L po 



2)r(^in|)re(i 



2L p a ) 



wpoL r(^ln2£ + l)r 6 (^ln^ 

\ L pa ' ) \ L p a 

distribution function decreases power-like: 



exp 



I- , 2p 
— — (In — 

2L V po 



L 



(43) 



For small values of p, when p < ^-exp(p-), the main contribution comes from the pole at the origin and the 



p 2 



l- 



2ttL 



(44) 



Thus the distribution function for the zeroth harmonic or (^-independent component has asymmetric form. 



B. Amplitude of the first harmonic contribution to the DOS moments 

By a Laplace transform with respect to x, Eq.(|3l|) with n = 1 is converted to 

(A + Sl m)Ql{m- A) - S m , = (to + l) 2 [Q\{m + 1; A) - Q\{m] A)] + to 2 [QJ(to - 1; A) - Q\{m- A)] (45) 



The 5 symbol on the left-hand side of this equation comes from the boundary condition (|32|). Eq. 
to the equation for the central part in the Berezinskii technique for strictly ID systems with open 
?7i 3> 1, this equation is transformed into a differential equation, 



corresponds 
30undaryH3 For 



i2 d 2 Qi (to; A) + 2m dQl(m ; A) 
dm 2 dm 



2to ^ ' ' - (A + s 1 m)Q* 1 (to; A) = (46) 



A change of the function to -&(z, A) = Q*(m; A), where z 1 = Amsi, reduces Eq.(E6|) to the Bessel equation 



d 2 $ | ld$ - | 1 + 4A ^_ Q ( 47 ) 
dz 2 z dz z 2 

Therefore Q\(m; A) can be expressed as 

Q\{m- A) - C 1 1 K 1+2g (2(ma 1 )^) (48) 
2(ms\) 2 

where q = - \ + ^J\+j. 

Eq. ([l8|) contains an unknown parameter C due to the neglect of the Kronecker symbol in Eq. (|4^) . 
On the other hand Eq.(^5|) has been solved by Mel'nikov^l who obtained the asymptotic solution of Q\(m\ A) for 
1 <C to <C s^ 1 as 

The comparison or Eq.([48|) with the asymptotic form (^9|) allows to determine C: 

C -4s« +1 r3(g + 1) (50) 

6 ~ 4Sl r(2g + 2)r( 2(Z + i) (50) 
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Taking the inverse Laplace transform, one obtains for Q*(m; x) 

/°° a\ l- a ( ^ 

^ e -|(A 2 +l)^ s ^ (mgl) -l K _. A(2(mSl) | ) r ^ 2 ^ (51) 

To get an expression for the first harmonic, (p k (e, L)) n= i, we substitute the solution ( |5l| ) into Eq. (|34|) , and sum over 
m, which can be done after the transformation of the sum into an integral over n — msi: 



</(e,Z)) n=1 = ^- > n \ > Sl k e «- / — Sl 2 e ' - 



2' kT^{k) °i i_ 00 2 7 r < ' 1 I*(-*A) 



dv _^ y2 r3(i^; 



2^ Sl 6 " r*(-iA') 



r(fc + iX + iX '- 1 )r{k + lX - lX '- 1 * 



(52) 



For convenience, we substitute below the variables A and A' by z and z' according to A = z + z' and X' — z — z' . 
The values of the integrals in Eq.(|52]) are determined by saddle points and poles. For 7 = ^~ In — < k — A, the 
contribution from the saddle point dominates: 



For 7 > fc — i the main contribution is given by the pole at z = i(k — A) and one gets 



(53) 



, k( nN ,p , fc 2(fc-i)r(2fc)r(2fc-i)r 6 (^) / ?- ^ (2fc _ lr _^. 



In contrast to the expression of (p (e, L)) n=0 for small dissipation, 77 — > [see Eq.(|42j)], the expression for (p (e, L)) n = 



increases strongly with si = 4t~?7 — > 0. It is illustrative to rewrite the prefactor of Eq. (p4j) as (2-7TpSi)2 — (87r 2 -^)2 
in terms of the level distance A = — j and the dissipation energy 77, the latter blurring the quantized energy levels. 
By decreasing 77, the energy levels are sharpened and the distribution function becomes a ^-function. 
Substituting Eq.(Q) into (|39|), one receives a normal logarithmic distribution for p > ^-exp(— p): 

^ ^ = Tl V ^ r(| + f m gfft m g - |)r^E ln g + ln § ex P L" 2Z ^ * F > \ (55) 

The logarithmic normal distribution function for the first harmonic is valid for a large range of p. Compar- 
ing Eq.(p5[) with Eq.(|43]) for winding number zero, it can be seen that Eq.(|55"|) contains in addition a prefactor 

= (itrjl^ p)~ x / 2 which increases with decreasing temperature. Thus the first harmonic increases with decreas- 



rir-p 

ing temperature faster in amplitude than the zeroth harmonic 



VI. CONCLUSION 



The distribution function for the local DOS in a one-channel ring threaded by a magnetic flux through the opening 
was studied in this paper. For this purpose, we constructed a new diagrammatic method as an extension of the 
Berezinskii techniqueEJ to the problem with periodic boundary conditions and in the presence of an external magnetic 



field. The equations obtained (fllCj) to ( |12| ) and (21) to ( J23| ) for the DOS and its fcth moments, respectively) are exact 
in the framework of the weak disorder limit kpl 3> 1. Eqs.([ll]) and ( |l2] ) are solved exactly, which gives the oscillation 
of p with the full flux for both weak localization and ballistic regimes. 

In contrast to the DOS problem, the equation for (p k (e, cf>; L)) is rather complicated and we succeeded to solve it 
for the weak localization limit when L > l^. In this limit, the leading contributions to arbitrary moments of the 
DOS oscillate with the halved period 4p The distribution functions for zeroth (insensitive to the magnetic field) 

and first (with a period of harmonics are calculated and logarithmic normal distributions (Eqs.(43) and (p 



are obtained for them, indicating large contributions from high moments of the DOS. For the zeroth harmonic, this 
normal logarithmic shape appears for the tail of the distribution, but for the first harmonic it covers the large range 
of p > 4^ exp(— p), i.e. the high moments give essential contributions not only on the tail but also in the vicinity of 
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the average value of the DOS. The distribution function for the first harmonic increases with decreasing the width of 
the energy levels or the dissipation parameter 77 (see Eq.(|55|)), which was introduced phenomenologically in the theory 
(Eq.(P). For r) — > 0, the distribution function W n= \(p) becomes a ^-function due to the quantization of the energy 
levels in the rings. The results for the DOS show that the amplitudes of all harmonics of p(e, <fi) are exponentially 
small in the weak localization regime [Eq. ([l4|)] , while the amplitudes of the higher moments in this regime [Eqs.([38|), 
(fll]), (42), and (jHJ)] are relatively large. Although we could not calculate higher moments of the DOS in the ballistic, 
regime, the amplitude of the average value of the DOS is large and seems to be consistent with experimental data.E-3 

It is also well known that the DOS of ID disordered crystalline systems is very sensitive to the filling factor. There 
exists disorder induced enhancement of the DOS for commensurable values of the electron wavelength A and the 
lattice constant a, when the electronic energy e satisfies the condition p(e) = k — ±1, ±2, • • • ± n and n — 2,3, ... , 
and the effect is pronounced for half-filling which corresponds to n = 2. The singularity in the DOS of ID disordered 
crystalline systems-jiear the middle of the band is known as a Dyson singularityE-3 which was studied for many ID 
electronic models.ElLa Notice that the Berezinskii method has also been applied to study the conductivity and the 
localization lengthEiea apartijrom the Dyson singularity in the middle of the band of a ID infinite lattice with both 
weak and strong disorder I 60 " 61 ! Our preliminary study shows that the real space diagrammatic method presented in 
this paper is applicable to study the Dyson singularity in the DOS of a ring for a half-filled energy band. This leads 
to a remarkable high amplitude of the persistent current as it is observed in the experiments, provided that the Peierls 
transition is suppressed by impurities and by weak transvers tunneling between the channels. 

The authors thank M. Kiselev for discussion. This work was supported by the SFB410. 



APPENDIX A: SOLUTION FOR Q (M,N;X) 



The first term on the right hand side of Eq. ( |ll|) can be removed through the transformation 
Q Q (m,n;x) = exp{-^ r (2m + n) 2 - jr m(m + n)}Qo(m, n ! x ) 



which gives for Eqs.(|ll|) and (|12j ) the simpler form 

j dQn (tTI, 71, x) o f ^ 1 -1 r*\ • 1 n \ f- 

l ; = — to expj — jQ (m — l,n + 2;x) with Q (m, n; x = V) = o m ,o 



d.r 



Laplace transformation from a; to A yields 



1 



\Qo(m, n; A) - (5 m ,o = -— m Qo(m — 1, n + 2;x — — ) 



Eq.([A3|) can be solved by iteration: 

Qo(m,n;X) 



l- 



V) 2 n 



/iA-jMj+n-l) 



where z 1 = Xl~ 



_ (ml) 2 r(z + i + ^)r(^ + i-z) 
~ a r(^ + TO + i + z)r(^ + TO + i-z) 

n ~ ' . The inverse Laplace transform gives for Qq 



r Q (m,n;x) = (to!) 2 ^ ex p{ + n ~ 1)}- 



3=0 



(-1) J U + n-2)\ 
j\(m — j)\ (m + j + n — 1)! 



(2j + n-l) 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



which, in connecti on w ith (Al), gives the final result Eq.([13|), where x is replaced by the full circumference L. The 
compliance of Eq.( A5) with the boundary condition is easily checked. Also, for to = 0, Qq(0, n; x) ~ 1, and for n = 
we get from the inverse Laplace transform of Eq.(A4) or from taking the limit of Eq. (|A5|) 



Q* (m, 0;x) = (1 - m - m£) + £ exp{^i(j + n-l)}( m ) (-1)^ Tn! g + n 2 | ! (2j + n - 1) (A6) 

3 — * 
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APPENDIX B: CALCULATION OF THE MIXING COEFFICIENT 



Using the relations 



dz - fe — 1 and / "' 



\z\<l 2ni 

we can transform Eq.(20) for <fii(m, m'; n + , n — n + ) to 



dz 



;1 2m z k+1 (l - z ) m - k+1 



(Bl) 



ifi(m, m; n , n — n ) 



dzi 1 / dz 2 1 / dz 3 1 / dz4 1 r (1 + Zi)(l + 22) — z 3 z 4 1 ' 



2tti z^ +1 / 2vri z ™' +1 / 2ni 7 " + + 1 T 2ni 



1 ( (L + z 1 )(L + z 2 )-z 3 z i y 
"- n+ +n(l-z 3 )(l-z 4 )-z 1 z 2 i [ ' 



Substituting z = Z1Z2, the dominant contribution for to 3> 1 comes from the pole at z = (1 — 23) (1 — Z4). Integrating 
over this new variable gives 



<pi(m, to ; n , n — n ) = 



(m + Z- 1)! / dz : 



1 



dza 1 



dZ4 



1 



(Z - l)!m! J 2m z m'-m+l+i J 2ni J 2m z 4 ~ n++1 

[(z 2 + (1 - z 3 )(l - z 4 ))(l + z 2 ) - z 2 z 3 z 4 ] 1 



(B3) 



(1 - z 3 ) m + l (l - z 4 ) m + l 
The remaining integrals are done in a similar way, resulting in 



(fi(m,m ;?i + ,rt - n + ) = - — - 



(2/)!(m + 1 + n+ - l)!(m + Z + n - n+ - 1)! 



(m' - m + Z)!(m - to' + Z)!(m + Z - l)!ra+!(n - n+)\(l - l)!m! 



(B4) 



Now we can collect all n + and n+ dependent terms in Eq. (pl|) and sum over n + and n+, introducing the mixing 
function <£>£ from Eq.(E8L). For large to, we can use Stirlings formula 



(m + o)! f,_ a , 

lim TTT m =1 

m->oo (m + fr)! 



(B5) 



to obtain 



$ fc (TO, to, to', to', n, n) = — - 



1 + e *o ) /fc 



(Z - l)!(fc - Z - l)!rt!n! \l J \m - m' + IJ \m -m' + I 
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2k -21 



(B6) 



Taking into account Q*(to,to, n;x) = Q*(m,n;x)d m ^n, this gives Eq. (|33|) . For to, to' 3> 1, the summations over Z 
and to' can be done. Following the paper [ plf|, we denote Ato = to — to', with Ato <C to for large to. The significant 
contributions to Eq.(p3[) come from Q*{m — Ato, n; a;) rs Q*(to, ra; a;). Hence we can rewrite Eq.(33) as 



(p (e,0;i)) =(— ) ^^[cos(2^— )J e =PJ — Q (to, n; x) 

m=On=0 ™° ' 



2/ 



2Zc - 2Z 



J J ^ \l + AmJ \k-l + AmJ (I- l)!(fc- I- 1)! 



1 



The last two sums result in 



2{k- 1) T 2 (2/c) 
fc(2fc- 1) T 5 (fc) 

which is used in the final result for the fcth moment of the DOS, Eq. 



(B7) 



(B8) 
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